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Abstract — The Welch-Berlekamp approach for Reed-Solomon 
(RS) codes forms a bridge between classical syndrome-based 
decoding algorithms and interpolation-based list-decoding pro- 
cedures for list size I = 1. It returns the univariate error-locator 
polynomial and the evaluation polynomial of the RS code as a 
j/-root. 

In this paper, we show the connection between the Welch- 
Berlekamp approach for a specific Interleaved Reed-Solomon 
code scheme and the Guruswami-Sudan principle. It turns out 
that the decoding of Interleaved RS codes can be formulated as a 
modified Guruswami-Sudan problem with a specific multiplicity 
assignment. We show that our new approach results in the same 
solution space as the Welch-Berlekamp scheme. Furthermore, 
we prove some important properties. 

Index Terms — Guruswami-Sudan (GS) interpolation, Reed- 
Solomon (RS) codes, Interleaved Reed-Solomon (IRS) codes 

I. Introduction 

The Guruswami-Sudan (GS) J6) approach for Reed- 
Solomon (RS) codes consists of an interpolation and a fac- 
torization step of a degree-restricted bivariate polynomial. 
The usage of multiplicities in the first stage improved the 
error-correcting capability of Sudan's original work |14|. 
The set of y-roots of the bivariate interpolation polynomial 
gives the candidates of the evaluation polynomials of the 
corresponding RS codes. The GS principle coincides with 
the Welch-Berlekamp (WB) approach [2| when the list size 
is £ = 1. Then, To = [(n — k)/2\ errors can be uniquely 
corrected, where n is the length and k the dimension of the 
RS code. 

Interleaved Reed-Solomon (IRS) codes are most effective if 
correlated errors affect all words of the interleaved scheme 
simultaneously (see H). Because of this, IRS codes are 
mainly considered in applications where error bursts occur. 
Bleichenbacher et al. Q, H formulated an IRS decoding 
procedure with the WB method. 

Our contribution covers the reformulation of the Bleichen- 
bacher approach in terms of a modified GS interpolation 
problem for a heterogeneous IRS scheme as it was investigated 
in |fl3l . The heterogeneous IRS code is built by virtual 
extension of an RS code. The rate restriction and the decoding 
radius of this scheme are comparable with the parameters of 
Sudan's original algorithm (where the multiplicity for each 
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point equals one). Also, the corresponding syndrome formu- 
lation (for Sudan done in fl2l . ifTTl ) is equivalent. Hence, it 
seems to be surprising that this scheme can be formulated as 
a modified GS interpolation problem, where the multiplicities 
are assigned in a specific manner. 

The paper is organized as follows. First, we shortly describe 



the GS principle for RS codes in Section III and outline 



important properties that we will use later. The connection to 
the WB approach is investigated in Section [TV] The virtual 
extension to an IRS code |fl3ll is described in Section [V] 
Section VI links the GS list-decoding procedure with the 
WB formulation of the previously described IRS scheme. 
Furthermore, the equivalence of both approaches is proved and 
an informal description is given. Finally, Section [VTT] concludes 
the paper. An example is given in the appendix. 

II. Definition and Notation 

Here and later, [n] denotes the set of integers {1, . . . , n} 
and [n]o denotes the set of integers {0, ...,n}. The entries 
of an to x n matrix S =|| Si.j | are denoted Sij, where 
i G [m — l]o and j G [n — l]o- A univariate polynomial of 
degree n is noted in the form A(x) = J^^o^iX 1 - A vector 
of length n is denoted by r = (ri,r2, . . . , r n ) T . 
Let at, Q!2) • • • j OL n be nonzero distinct elements (code lo- 
cators) of the finite field F = GF(q) of size q. C = 
{ai, . . • , a n } is the set containing all code locators. Denote 

f(C) = (f(a i ),...J(a n )) 

for a given polynomial f(x) over F . 

An RS code lZS(n, k) over F with n < q is given by 

1lS{n,k)={c=f{C) :f(x)e¥ k [x}}, (1) 

where F^[x] stands for the set of all univariate polynomials 
with degree less than k and indeterminate x. 
RS codes are known to be maximum distance separable 
(MDS), i.e., their minimum Hamming distance is d = n—k+1. 

III. The GS principle and the univariate 
formulation 

A. Guruswami-Sudan Approach for Reed-Solomon Codes 

Let the n points {(ai, rj)}™ =1 , where ai, r.i G F and r = 
(r 1? . . . ,r n ) denotes the received word, be interpolated by a 
bivariate polynomial Q(x,y). The number of errors, that can 
be corrected, is denoted by r. The parameter s is the order of 



Dd 



(2) 



multiplicity of the bivariate interpolation polynomial in the 
GS algorithm. The list size is denoted by I. The nonzero 
interpolation polynomial Q(x, y) has to satisfy the following 
degree conditions: 

deg 01 Q(x,y) < £, 
_deg lk _ 1 Q(x,y) < s(n - r) 

where deg u v a(x, y) — ud x + vd y is the (it, w)-weighted de- 
gree of a bivariate polynomial a(x, y) = J2i=o 2j=o a i,j xl y J ■ 
The interpolation constraints are: 

Id := [Q [a ' b] (a^ n) = Vi G [n] andVa + b < s] , (3) 

where Q^ a ' b \x,y) represents the mixed Hasse derivative 
(see [7.| for definition) of the polynomial Q(x,y) G F[x, j/] . 
Analogously, one can say that the GS polynomial must have 
a multiplicity of s at each point (dfj, rj). 

B. Univariate Formulation of Gurusw ami-Sudan 

In (TJ, ifTTl the univariate reformulation of the bivariate GS 
interpolation problem (£ey equations) was derived. Here, we 
state some basic properties that will be used later on. 

Proposition 1 (Augot-Zeh |1|) Given s > 1, let Q(x,y) = 
Y] t —Q {x)y t be the Gurusw ami- Sudan interpolation poly- 
nomial that satisfies (|2j and (jij and let R(x) be the Lagrange 
interpolation polynomial, such that R{ai) — Ti Vi £ [n] 
holds. Furthermore, let G{x) = Y\^—^{x — ctj). Then, Q(x,y) 
satisfies (jij, if and only if there exist s polynomials B^ b \x) £ 
¥[x] V6 G [s - 1] with: 



Q [b] (x,R(x))=B^(x)-G(x) s - b , 
where deg (x) < l(n — k) — st + b. 



(4) 



We remark that Q^{x,y) := Q^°^(x,y) denotes the 6-th 
Hasse derivative of the bivariate polynomial Q(x,y) with 
respect to the variable y. 

IV. Welch-Berlekamp approach as List-1 Decoder 

We recall a simplified version (as in J5) or J8] Ch. 5]) of 
the WB approach iflOl Ch. 7.2] (2| for decoding RS codes up 
to half the minimum distance (tq = [(n — k)/2\). It is seen 
as special case of the list-decoding problem of GS. 
The interpolation polynomial Q(x,y) of the GS algorithm for 
£ — s — 1 has the following form: 

Q(x,y) = Q^(x) + Q^(x)y, 

where deg (x) < n — r and deg (x) < n — T — k + l. 
Condition ^ simplifies to Q(oei, r{) — V« G [n] and gives n 
linear equations. The codeword c coincides with the received 
word r in at least n — r positions. Therefore, we have: 

Q(x,.f(x)) = Q^(x) + f(x)-QW(x) = 0. 

So f(x) — ~Q^(x)/Q^(x) and we can rewrite the original 
interpolation polynomial: 

Q^(xY 



Q(x,y) 



)(!) 



(x)-[y + 



QW(x) 



Q^(x)-(y-f(x)). 



Clearly, QW(i) is the error-locator polynomial (ELP), be- 
cause it vanishes for tq a/s. Let the classical ELP A(a;) = 
rijej r ( :E — a j)> where J is the set of error locations. Then, 
we can write: 



Q{x,y) = A(x) ■ (y - f(x)). 



(5) 



In the WB decoding procedure the polynomial Q(x, y) of |5]) 
is determined by solving n linear homogeneous equations. The 
standard syndrome-based decoding procedure, that consists of 
To equations for the ELP, can be derived by reducing the WB 
equation. 

V. Virtual Extension to an IRS code 
A. Basic Principle 

We shortly describe the Schmidt-Sidorenko-Bossert 
scheme |[T3l where an RS code is virtually extended to an 
IRS code. This IRS code is denoted by VllZS(n, k, s), 
where n and k are the original parameters of the 7ZS(n, k) 
code. The parameter s denotes the order of interleaving. Let 
P( x ) = Y^j=o PjX'' he a univariate polynomial in F„[a;]. 
Then, 

n-l 



P 



'(X) 



3=0 



is the polynomial in ¥ n [x] where each coefficient is raised 
to the power i. Analogously, c <l> denotes the vector 
(c\, . . . , c*J T . The virtual IRS code can be defined as follows. 



Definition 1 (Virtual Extension to an IRS code |13|) Let 

lZS(n, k) be an RS code with the evaluation polynomials 
f(x) as defined in ([7J. The virtually extended Interleaved 
Reed-Solomon code VTlZS(n, k, s) of order s is given by 



V2US{n, k, s) = 



/C<!>\ 
,<2> 



//(£) 



:f(x)e¥ k [x] \ 
■■ (/(*)) 2 eF 3(t _ 1)+1 [i] 



\f s (C) :(f(x)y e¥ s(k _ 1)+1 [x]J 

Clearly, the parameter s must satisfy s(k — 1) + 1 < n. The 
scheme is restricted to low-rate RS codes and allows to decode 
beyond half the minimum distance. The virtual extension is 
illustrated in Figure [T] where the information length of the 
i-th codeword is fcW — i(k — 1) + 1. The decoding procedure 
for the virtual extension of an RS code is as follows; the 
elements of received word r = c + e are raised to the power 



. , s (r 



<2> r <3> 



, r <s> ) and a heterogeneous 



IRS code is obtained. Clearly, through the virtual extension, 
the error is also "extended" and every single received word 
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Fig. 1. Illustration of an 7?.cS(31,4) code that has been virtually extended 
with interleaving factor s = 3. The errors in the 7?.cS(31, 4) code are extended 
to burst errors in the VI7?.cS(31, 4, 3) code. 



v <l> is erroneous at the same positions. Due to the additional 
equations, the decoding radius is increased to: 



(s+l\ 
\ 2 



1) 



1 



(6) 



The radius r is greater than To = [(n — k)/2\ for RS codes 
with code rate R < 1/3. (For further details (e.g. increased 
failure probability) of this scheme, see (TSJ). We remark that 
the rate-restriction and the increased decoding radius coincide 
with the original Sudan algorithm (where the multiplicity s 
equals one for all points (a^r;)). Nevertheless, we will show 
that this scheme is equivalent to a GS interpolation problem 
with a modified multiplicity assignment and stricter degree 
constraints. To start the logical chain, we will describe in the 
following the corresponding system of equations of the s WB 
equations for a VllZS(n, k, s) code. 

B. Matrix form of the Set of Equations 

Bleichenbacher et al. I3, (4) described the WB formulation 
for IRS codes. We recall this approach for the virtually 
extended Reed-Solomon code VllZS(n, k, s). 
Clearly, we have s WB-equations (see |5]0 of the form: 

Q< b >(x,y) = A(x)-(y b ~f»(x)) 

=:Q^(x)y b -Q^(x), (7) 

for all b £ [s - 1] . 

For every single WB polynomial Q <b> (cti,ri) = holds 
(i £ [n]). Note, that through the virtual extension, each 
received word r <z> has its errors at the same position and 
therefore we search one common ELP A(x). We represent the 
sn constraints of system |7]i in matrix form. Therefore, let the 
n x (r + i(k — 1) + 1) matrix Mi be: 



/I 


Ql 


al ■ 


Ni-l 
■ "l 


1 


a- 2 




Ni-1 

■ a 2 


1 


<*3 


°§ • 


„JV<-1 
• a 3 



or. 



(8) 



where JVj 

N = 



i=0 



N; 



i ■ (k — 1) + 1 and let N be defined as: 

= (s + l)(r + l)+ ( S+ ^)(k-l). (9) 



Furthermore, let the nxn matrix R have the following form: 
(r x ■■■ 



R 



r 2 



\0 







(10) 



Now, we can write the s polynomial equations from (j7]) in 
matrix notation. Let Q = (Q (0 \ Q (1 \ . . . , Q (s:) ) T , where 
Q(0 = (Q^ ) ) Qi i) ,...,QS (5 _, )(fe _ 1) ) T . The homogeneous 
set of equations is of the form A • Q = 0, where the sn x N 
matrix A is: 



( ° 




A = 





-M 2 



-Mi R- Mo \ 
R 2 M 



(11) 



\-M s • • • R s • M / 
The vector gives the coefficients of the ELP A(x). 
VI. Reformulation as a modified 

GURUSWAMI-SUDAN PROBLEM 

A. Specific Multiplicity Assignment 

In this section, we formulate the decoding of an lZS(n, k) 
code virtually extended to a V17ZS(n, k, s) code as a modified 
GS interpolation problem. The constraints of the bivariate 
interpolation polynomial with multiplicities are a modified 
version of the general GS algorithm introduced in Section III 
We show the corresponding homogeneous set of equations 
and prove the equivalence to the one of Bleichenbacher et 
al. (see (fTl)). 

Let Q(x, y) be a bivariate polynomial of ¥[x, y] \{0}, where 



DCo := 



deg 0jl Q{x,y) < s, 



?(*)/ 



degQ w O) <T+{s-t)-(k-l) 
The modified interpolation constraints for Q(x, y) are: 
IC 2 



(12) 



Q l 1 ( ai , n)=0 Vi£ [n] and V6 £ [s - 1] , 

(13) 

where the parameter s is such that s(k — 1) + 1 < n holds 
and Q^ b \x,y) denotes the 6-th Hasse derivative with respect 
to the variable y of the polynomial Q(x,y). 

Theorem 1 // exists at least one nonzero polynomial Q(x, y) 



which satisfies conditions (13 'I. 



Proof: Condition ( p"3] > gives sn homogeneous linear equa- 
tions to the coefficients. The number of possible coefficients 
is TV (as defined in (|9|), therefore we get a nonzero solution 
for the decoding radius r as in Equation |6]) 



The bivariate polynomial Q(x,y) that fulfills condition (13 1 
has multiplicity s for all n — r error-free positions and 



multiplicity one for all r error positions. Let us state this 
property in the following theorem. 

Theorem 2 The bivariate polynomial Q(x, y) under the con- 
straints DC% and IC2 can be written as: 



Q(x,y) = Q {s \x)-(y-f(x)y, 



(14) 



where Q (x) is the ELP and f(x) is the information poly- 
nomial of the RS code (see definition ([TJj. 



Proof: Let us consider the "last" (s — l)-th Hasse 
derivative of Q(x, y) with respect to the variable y: 



Q (s \x)y 



— Q s 1 (a) + s ■ Q s (x)y 

which is by ( fl3| ) zero for the set {(c^, ri)}™ =1 . Clearly, 
g} S 1 '(.T,y) is a WB polynomial for the lZS(n, k) code with 

(s— 1) (s) 

information polynomial fix) = —Q (x)/s ■ Q (x) (see 
Section HVl 



The (s — 2)-th Hasse derivative of the interpolation polynomial 
q' S 2 '(x, y) can now be rewritten as: 



0,2/) 



s 

s-2 



(*) + 



s- 1 
s-2 



=Q is 2 \x) + ( s -l)-Q {s l \x)y- 



S (s-1)-Q (s \x)y 2 
s(s-l).Q {s \x)-(y 2 -f(x)y)+ 



Q [x) 



(15) 



where 



\x,f(x)) = 



from the interpolation constraints holds. We can now express 

Q {s ~ 2 \x) as: 



(x) 



\ S {s-l)-Q {s \x)-{f(xf-2f{xf) 



l -s(s-l)-f{xf-Q {s \x). 



(16) 



Substituting this into ( 15 1, we obtain for the (s — 2)-th Hasse 
derivative of Q(x,y): 



*[*-2] 



1 



(x,y) = ~s(s-l)-Q^(x)-(y-f(x)) 



which has multiplicity two at the n — r error-free positions 
and multiplicity one at the r erroneous positions. By induction 
we can state that (y — f(x)) s \Q(x,y) and Q^ s \x)\Q(x,y). 
From DC2 we know, that no other polynomial factor occurs 
in Q(x,y). ■ 



B. Informal Description 

The degree condition DC2 and the interpolation constraint 
IC2 for the polynomial Q(x, y) are a subset of the general 
GS list-decoding constraints DC\ and IC\. The y-degree 
of Q(x,y) corresponds to the number of codewords of the 
VIlZS(n, k, s) code. Similar to the univariate formulation of 
the original GS interpolation problem (see Q) it is sufficient 
to consider only the Hasse derivatives with respect to variable 

y- 

In the original GS algorithm the 6-th Hasse derivative of the 
interpolation polynomial Q(x,y) is divisible by G(x) l - S ^ b \ 
where G(x) = Y[i=i( x ~ a i) an( l 71 denotes the code 
length. In our case the 6-th Hasse derivative of the modified 
interpolation polynomial Q(x,y) is divisible by G(x) l - S ^ b \ 

Here, G(x) — riie[n]\j r ( :c — a i) an< ^ W \ 3 ^ s trie se * °^ 
error-free positions. 

( s \ t s \ 

Furthermore, the ELP Q (x), where degQ (x) can be 

greater than [On — k)/2\, is a factor of Q(x,y). The zeros 
of Q (x) have multiplicity one. 

The scheme of Section [V] virtually extends the received vector 
r = (r 1 ,r 2 , ■ ■ ■ ,r n ) of an 7ZS(n,k) code to s received 
words r <%> — (r\, r|, • • ■ , r^) Vi € [s] of s different 
lZS(n,i(k — 1) + 1) codes with equal code length n. 

C. Set of Equations 



Now, we consider the homogeneous set of equations (13 1. 



We have B • Q = 0, where Q is the vector notation_of the 
interpolation polynomial Q(x,y). The sn x N matrix B can 
be written as: 



/ 


V M s RM S _! 







;jRMi 



R s_1 Mi 



gRMo \ 
,)R 2 M 



R S M / 



where the sub-matrices Mj and R are defined in (|8]l and ( fTO) . 
The binomial coefficients come from the Hasse derivatives of 

Q(x,y): 



7#1 



(x,y) = 



t=b 



%Q (t \x)y^. 



(17) 



Note, that the first n rows of matrix B correspond to the 
(s — l)th Hasse derivative of the polynomial Q(x,y). The 
second n rows represents the n interpolation constraints of 
the (s — 2)-th Hasse derivative and so on. In the last n rows 
of matrix B the interpolation polynomial Q(x, y) occurs with 
all terms. 

D. Equivalence of Both Sets of Equations 

In the following, we show the equivalence between the 
systems of equations determining the IRS scheme of Section[V] 
and the one determining the modified GS interpolation poly- 
nomial Q(x,y). Due to space limitations we will sketch the 
basic steps of the proof. 



First, let us consider the relation between vectors Q and Q: 

Q(x,y) =Q (s \x)(y-f(x)y 



si=0 



=Q w (*)-(£( :)(-i)y-7(*) 

In vector notation, we have: 

(Q (0) ,.. 



Q (s - 2 \ 
(o) 



s — 1, 

(s-2) ^r(s-l) ^(ir 



s -2 

(Q W ,...,Q V " ",Q V " ".Q"") 
Let the matrix B be such that: 

B Q = B Q. 
Matrix B is then (first column not printed): 



B = 



)M 2 



,)R 2 M ' 



V C-0 R ' M I?M, / 

After simplification, we obtain: 

B = 

/ ••• -Mi RM \ 

|s(s-l)M 2 —s(s — 1)RM] |s(s - 1)R 2 M 



\ ••• ••• -sR s_1 Mi R S M / 

The second band of n rows of matrix B can be multiplied 
with — Rs(s — l)-times the first band of B and then divided 
by — |s(s— 1). We obtain the second band of n rows of matrix 
A. Repeating this operation, matrix B can be transformed into 
matrix A (jTTJ. 

VII. Conclusion 

We investigated a virtual extension of an RS code to an 
IRS code from an interpolation-based list-decoding approach 
point of view. 

The Bleichenbacher scheme was used to form the system 
of equations for the IRS scheme (based on a virtual exten- 
sion). Then, the original constraints of the GS list-decoding 
algorithm were modified and the equivalence of the resulting 
system of equations with the Bleichenbacher scheme for the 
IRS code has been shown. 
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Appendix 

Let us consider an 7?.iS(16,4) code over F = GF (17) 
with parameter s — 2 (number of interleaving and multiplicity 
for the modified GS algorithm). The corresponding increased 
decoding radius is r — 7 (see d6])). 

The code locators are ai = a Vz € [n], where a is 3. For 
the information polynomial f(x) = 1 + x + x 2 + x 3 (see (Jl]i) 
and an error e of weight r = 7 we get the following vectors: 

c = (4, 6, 4, 6, 0, 3, 12, 2, 0, 14, 7, 9, 0, 15, 15, 4) 
e = (1, 2, 3, 4, 5, 6, 7, 0, 0, 0, 0, 0, 0, 0, 0, 0) 
r <1> = (5, 8, 7, 10, 5, 9, 2, 2, 0, 14, 7, 9, 0, 15, 15, 4) 



„<2> 



(8, 13, 15, 15, 8, 13, 4, 4, 0, 9, 15, 13, 0, 4, 4, 16). 



The conditions ( p~3] > on the modified bivariate polynomial 
Q(x,y) give the following solution: 

Q =(5, 14, 8, 6, 14, 9, 5, 9, 12, 12, 4, 2, 3, 16, 5, 9, 11, 
15, 13, 4, 7, 2, 16, 4, 16, 5, 4, 2, 4, 3, 12, 5, 16). 

And the corresponding modified bivariate interpolation poly- 
nomial; 

Q(x, y) ={x + 2)(x + 4) (as + 7){x + 8)(x + 12)0 + 14)- 
(x + 16) (y + I6x 3 + 16x 2 + 16x + 16) 2 , 

is factorizable as stated in Theorem [2] 



